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Abstract 


An  automated  method  has  been  implemented  in  MATLAB®  to 
compensate  for  signal  dispersion  in  antenna  structures.  We  have 
explored  postprocessing  techniques  that  involve  frequency 
transforms  and  deconvolution.  The  method  has  been  applied  to 
transient  signals  measured  from  a  variety  of  different  antennas  and 
impulse  sources.  The  technique  has  proved  to  be  a  valuable  tool  in 
reconstructing  fast  transient  signals  with  inexpensive  high-gain 
log-periodic  antennas  instead  of  more  expensive,  high  fidelity 
wideband  horns. 
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1.  Background 


High  fidelity  wideband  antennas  have  been  developed  to  preserve  the  pulse 
shape  of  ultra  wideband  transient  waveforms  [1].  These  antennas  are  typi¬ 
cally  horn  structures  that  occupy  large  volumes  with  low  gain.  The  gain  is 
lower  for  the  longer  wavelengths  because  the  gain  must  vary  linearly  with 
frequency  to  maintain  a  constant  effective  length.  These  anteimas  have  been 
very  effective  for  large  signal  strength,  high  fidelity  data-collection 
requirements. 

Commercially  available,  high-gain  log-periodic  (LP)  antennas  would  be 
useful  for  capturing  small  signal  transients  that  are  below  the  minimum 
detectable  signal  level  of  the  bulkier  large  aperture  antennas.  The  higher 
gain  LP  antennas  act  as  filters  that  modify  the  phase  and  amplitude  of  the 
input  signal  as  a  function  of  frequency.  Every  antenna  disperses  a  pulsed 
signal  to  some  extent.  Some  antermas  are  more  dispersive  than  others  be¬ 
cause  of  materials  used,  construction,  or  both.  The  purpose  of  this  work  is 
to  correct  for  the  dispersion  with  digital  filter  postprocessing  and  to  recon¬ 
struct  the  original  signal. 

Dispersion  from  an  antenna  is  characterized  by  group  delay  and  amplitude 
distortion.  Every  antenna  has  a  characteristic  dispersion.  Some  antennas 
are  large  physical  structures  that  delay  the  signal  because  they  have  many 
different  length-resonant  elements  (as  in  a  log-periodic  antenna);  others  in¬ 
troduce  signal  dispersion  because  cable  characteristics  and/or  baluns  are 
added  to  the  channel  path.  Figure  1  shows  the  dispersion  of  a  pulse  caused 
by  a  log-periodic  antenna. 


Figure  1.  (a)  The  direct  (a)  0,1 

output  of  a  Picosecond 
Pulse  Labs  4015C  0 

pulse  generator  and 
(b)  as  transmitted 
between  a  pair  of  log- 
periodic  antennas. 
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2.  Approach 


The  approach  is  to  measure  the  response  of  a  pair  (transmitting/ receiving) 
of  identical  antennas  to  a  wideband  source.  An  impulse  source  generates  a 
signal  that  is  transmitted  through  a  radiating  antenna  placed  at  one  end  of 
an  anechoic  chamber.  To  obtain  the  true  response  of  the  antenna,  a  very 
wideband  (0  to  10  GHz)  pulser  is  used.  A  receiving  antenna  is  cormected  to 
a  digital  oscilloscope  at  the  other  end  of  the  anechoic  chamber.  The  digital 
oscilloscope  acts  as  a  base-band  receiver.  Data  acquired  by  the  digital  oscil¬ 
loscope  are  downloaded  to  a  portable  computer,  where  they  are  filtered 
and  transformed  via  digital  signal  processing  (DSP)  in  MATLAB®  [2]  as 
described  next.  In  addition,  the  source  pulse  is  entered  directly  into  the  digital 
oscilloscope  and  recorded.  These  two  signals  are  deconvolved  to  obtain  the 
antenna  response  in  the  following  manner: 

_FFr(|w.) 

‘^TransFunc.FFT  cct/c  t  ’ 
hH^OThruj 

in  which 

^Direct  is  the  direct  voltage  signal  from  the  pulser, 

Sxhru  is  the  response  of  the  pulser  through  the  antennas  of  interest, 
SiransFuncjFT  is  the  frequency  response  of  the  transfer  function,  and 
FFT  is  the  fast  Fourier  transform. 

The  transfer  function  response  filter  is  calculated  from  the  ratio  of  the  direct 
and  transmitted  signal.  It  is  then  applied  to  the  data  taken  with  the  original 
or  other  pulsed  sources.  By  applying  the  transfer  fimction  to  the  measured 
data,  we  can  remove  the  dispersive  effect  of  the  antennas.  This  is  done  m 
the  following  manner: 


^  _  PP^(‘^Measured) 

(2) 

‘^CorrectedlFFT  ~  f  ^  ^ 

\^‘^TransFunc,FFTj 

Co^ttAFFT  =  lFn(sco,recte<lFFT) 

(3) 

in  which 

^Measured  is  the  measured  signal  from  the  receiving  antenna, 

Scorrected  is  the  Corrected  time  response  of  the  pulser  through  the  antermas 
of  interest, 

SjransFunc.FFT  is  the  frequency  response  of  the  transfer  function,  and 
Ifft  is  the  inverse  fast  Fourier  transform  (FFT). 
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3.  Experimental  Configuration 


Figure  2  shows  the  experimental  arrangement.  The  configuration  consists 
of  a  source  for  generating  fast  pulses  connected  to  a  transmitting  antenna 
and  a  Tektronix  11801C  sampling  oscilloscope  connected  to  a  receiving  an¬ 
tenna.  The  11801C  sampling  oscilloscope  uses  an  SD-24  sampling  head  that 
is  capable  of  measurements  as  great  as  20  GHz.  The  sampling  oscilloscope 
is  connected  to  the  source  via  a  cable  for  triggering. 

Several  impulse  sources  have  been  used  during  the  investigations  to  gener¬ 
ate  the  wideband  signals.  These  include  (1)  the  PicoSecond  Pulse  Labs  (PSPL) 
4015C  pulser,  which  has  a  nominal  rise  time  of  17ps  and  peak  output  of  9  V; 
(2)  the  Avtech  pulser,  which  supports  variable  rise  time  and  pulse  width; 
and  (3)  the  hand-held  (HH2)  generator,  which  is  a  reed  switch  pulser  built 
in  the  Amy  Research  Laboratory  for  impulse  applications  [3].  These  sources 
are  listed  in  table  1. 

Antenna  pairs  that  were  evaluated  and  compared  to  each  other  include  a 
commercial  log-periodic  antenna,  an  EMCO  double-ridge  wave  guide  horn, 
and  two  versions  of  a  resistively  loaded  transverse  electromagnetic  horn 
labeled  as  T2  and  T3.  The  characteristics  of  these  antermas  were  measured 
with  identical  pairs  of  antennas  during  these  experiments  and  are  listed  in 
table  2. 


Figure  2. 

Experimental 

schematic. 
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Table  1.  Sources 

Source 

Rise  time 

Output 

(V) 

Repetition 

rate 

PicoSecond  Pulse  Labs  4015C 

17  ps,  nominal 

9 

10  kHz 

Hand-held  pulser  (HH2) 

100  ps 

700 

100  Hz 

Avtech  pulser 

5  to  100  ns 

300 

IkHz 

Table  2.  Antennas 

Antenna 

Bandwidth 

Gain 

(dB) 

T3 

20  MHz  to  8  GHz 

-2 

T2 

100  MHz  to  4  GHz 

-5 

Log  Periodic 

100  MHz  to  1.3  GHz 

8 

EMCO 

200  MHz  to  2  GHz 

12 

5 


4.  Procedure 


The  transfer  function  technique  described  in  section  2  was  applied  to  two 
different  pulsed  sources  and  four  different  anterma  pairs.  As  mentioned 
earlier,  two  waveforms  are  required:  (1)  the  signal  directly  from  the  im¬ 
pulse  source,  and  (2)  the  signal  from  the  source  through  the  antenna  pair. 
The  signal's  input  to  and  output  from  the  T3  anterma  pair  (for  example)  are 
shown  in  figure  3.  The  derivative  of  the  "direct"  pulse  is  taken.  This  is  done 
because  an  anterma  transmits  according  to  the  derivative  of  the  signal  fed 
into  it.  Moving  charges  are  required  to  generate  the  crossed  fields  in  the 
anterma.  A  DC  signal  will  not  radiate.  Thus,  the  derivative  represents  the 
ideal  behavior  of  radiation  generated  in  an  anterma.  However,  since  no  an¬ 
terma  has  an  infinite  bandwidth,  signals  through  them  will  be  limited  in 
bandwidth,  as  shown  in  figure  3. 

The  FFT  is  calculated  for  the  two  waveforms  (direct  and  through  anteimas) 
shown  in  figure  4.  We  then  obtain  the  transfer  by  dividing  the  two  FFT 
signals,  as  described  in  equation  (1).  This  transfer  fimction  (see  fig.  5),  which 
is  in  the  frequency  domain,  is  then  plotted  and  saved. 

Figure  6  shows  the  complete  flow  chart  for  this  process.  This  flow  chart  has 
been  converted  to  a  MATLAB  [2]  script  file.  The  code  is  shown  in  appendix  A. 
Measurements  of  "direct"  and  "through"  are  saved.  These  files  are  input  to 
the  code  through  a  series  of  menus.  Any  attenuation  or  amplification  used 
during  the  measurements  is  also  input  to  obtain  proper  scaling  of  the  wave¬ 
forms.  After  the  user  enters  these  data,  the  transfer  fimction  is  cornputed 
and  saved  in  ASCII  format  for  later  use  in  the  correction  process. 

This  transfer  function  is  then  used  as  a  filter  and  applied  to  signals  received 
through  the  antermas  from  other  sources.  Specifically,  the  signal  transmit¬ 
ted  through  the  anterma  pair,  which  will  be  referred  to  as  the  "thru"  signal, 
is  read.  An  FFT  is  taken  of  the  "thru"  signal.  The  transfer  function  is  also 
read  and  interpolated,  if  necessary,  to  match  the  sample  interval  and  length 
of  the  "thru"  waveform.  The  two  signals  are  then  divided  as  described  in 
equation  (2)  to  create  a  corrected  frequency-domain  file.  This  file  is  then 
filtered  with  a  band-pass  filter  to  eliminate  high  and  low  frequency  noise 
and  cmomalies  that  are  outside  the  receiving  anterma's  operating  band.  These 
anomalies  are  spuriously  generated  during  the  compensation  process.  Then, 
an  inverse  fast  Fourier  transform  is  applied  to  obtain  the  corrected  time- 
domain  signal.  It  is  then  plotted  and  saved.  Figure  7  shows  the  flow  chart 
for  this  process.  The  code  that  performs  these  operations  is  also  shown  in 
appendix  A.  The  program  uses  a  series  of  menus  to  guide  the  user  through 
the  entire  process. 
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Figures.  Signals  input 
to  (green)  and  output 
from  (blue)  an  antenna 
pair. 


lesion. 


Figure  4.  Fast  Fourier 
transforms  (FFTs)  of 
the  input  and  output 
signals  of  figure  3. 
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Figures.  Transfer 
function  generated 
from  input  and  output 
waveforms. 


Figure  6.  Flowchart  of 
process  to  generate 
transfer  function. 
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Figure  7.  Flowchart  of 
signal  correction 
process. 


The  effectiveness  of  the  compensation  algorithm  varied  for  the  multiple 
antenna-source  combinations.  Figure  8  shows  relatively  good  reconstruc¬ 
tion  when  the  hand-held  source  and  T2  antenna  pair  were  used.  The  first 
plot  in  each  figure  shows  (1)  signal  through  the  antenna  pair  (blue),  (2)  com¬ 
pensated  version  of  this  signal  (green),  and  (3)  derivative  of  the  signal  taken 
directly  from  the  output  of  the  pulser  (red).  The  rise  times  of  these  signals 
are  included  in  the  legend  for  comparison.  The  second  plot  shows  the  fre¬ 
quency  responses  (FFTs)  of  these  signals.  Overall  effectiveness  was  based 
on  how  close  the  reconstructed  or  compensated  signal  was  to  the  derivative 
of  the  direct  signal.  Factors  used  to  compare  the  signals  included  the  rise 
times,  pulse  widths,  and  pulse  shapes  in  the  time  domain,  and  the  frequency 
content  and  amplitudes  in  the  frequency  domain.  As  the  sample  rates  (and 
Af)  vary  from  the  examples  shown,  so  does  the  calculated  derivation  of  the 
input  voltage.  This  variation  is  a  natural  result  of  finite  differences  inherent 
in  DSP. 

Figure  9  shows  the  reconstruction  when  the  hand-held  source  was  used 
with  the  EMCO  antermas.  Based  on  the  rise  times,  pulse  widths,  and  pulse 
shapes  of  the  compensated  and  original  signals,  very  good  signal  recon¬ 
struction  was  achieved  in  this  case. 

Figures  10  and  11  show  lower  quality  reconstruction  of  the  received  signal. 
One  of  the  main  reasons  for  the  poor  compensation  was  that  our  enclosed 
test  chamber  prevented  adequate  propagation  of  signals  below  200  MHz. 
At  200  MHz,  a  significant  ringing  was  apparent  in  some  of  the  signals,  which 
was  also  attributable  to  the  use  of  an  enclosed  metal  chamber.  We  believe 
this  is  why  the  results  were  poor  for  the  Avtech  pulser,  which  consists  pri¬ 
marily  of  signals  below  200  MHz.  Results  with  other  source-anterma  com¬ 
binations  are  shown  in  appendix  B. 

The  width  of  pulses  produced  by  the  Avtech  source  varied  from  pulse  to 
pulse.  As  a  result,  the  oscilloscope's  measuring  method  of  sampling  pro¬ 
duced  a  pulse  with  a  distorted  trailing  edge.  This  resulted  in  erratic  fre¬ 
quency  transforms  because  of  jitter  in  the  trailing  edge  of  the  pulse,  which 
is  shown  in  figure  12. 

The  compensatory  method  used  signals  that  were  collected  with  different 
sampling  rates  and  recording  lengths.  The  effectiveness  of  these  multirate 
compensations  varied  throughout  the  combinations  of  antennas,  sources, 
and  recording  devices.  The  complete  data  set  and  results  of  its  corrections 
are  shown  in  appendix  B. 
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Figure  8. 

Reconstruction  of 
signal  from  HH2 
source  through  T2 
antenna  pair  (see  text 
for  discussion). 
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Figure  9. 

Reconstruction  of 
signal  from  HH2 
source  through  EMCO 
antenna  pair. 
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Figure  10. 
Reconstruction  of 
signal  from  Avtech 
source  through  T2 
antenna  pair. 
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Figure  11. 
Reconstruction  of 
signal  from  Avtech 
source  through  EMCO 
antenna  pair. 
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Figure  12.  Jitter  in 
Avtech  source. 
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6.  Conclusions 


Wideband  sources  and  antennas  provide  better  compensation.  Care  must 
be  taken  to  ensure  that  the  entire  bandwidth  range  is  covered  when  one  is 
obtaining  the  antenna  response,  that  is,  the  reference  measurement,  espe¬ 
cially  for  the  lower  frequencies.  Our  instruments  limited  us  to  1,024  samples, 
which  may  not  have  been  adequate  to  calculate  the  low  frequency  response. 
With  the  latest  versions  of  fast  analog-to-ditigal  converters  and  megabyte 
record  lengths  available,  this  problem  should  be  eliminated. 

Future  work  being  planned  will  include  correcting  for  measurements  of 
antennas  that  are  off  boresight,  correcting  for  measurements  performed 
outside  the  normal  operating  frequency  ranges,  using  full  discrete  Fourier 
transforms  instead  of  FFTs  for  signal  correction,  and  using  other  digital  fil¬ 
tering  in  the  correction  process. 
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Appendix  A.  Results 


This  appendix  includes  the  internally  documented  script  files.  The  MATLAB 
[2]  code  developed  for  these  calculations  is  shown  below. 

%  comp_tranFuncGen.m 
%  4/99  N  Tesny 

%  generates  FFT  transfer  function  of  two  pulser  data  signals 
%  this  program  doesn't  handle  differing  sample  rates  yet 

FALSE=0 ;  TRUE=:~FALSE  ; 

%filpre  =  'g: /uwb/other_subjects/pspl/ ' 
filpre  =  'g:\uwb\other_subjects\lpDecon\'; 

%  ask  user  name  of  file  for  thru: 
pname=^filpre; 

[ f ilthru, pname]  =  uigetf ile ( [pname asc Select  THRU 
file' ) 

if  ( (filthru{l) -=0)  &  { filthru { length ( fi Ithru )) ~=0 )) , 
path_thru=pname ; 
else 

pa  th__thru= f  i  Ipre  ; 

end 

%  enter  attenuation  for  THRU  file: 

OKflag=FALSE; 
while  OKf lag“=FALSE , 

prompt  =  ('Enter  dB  attenuation Enter  ratio  attenua¬ 
tion:  '  }  ; 

titlel  =  'Input  Attenuation  for  THRU  file'; 
lines=  1; 

def  =  { '20' , '1' } ; 

answer  =  inputdlg (prompt , titlel , lines , def ) 

f =cell2struct (answer (1) , 'd' , 1) ; 
als=f . d; 

al=str2num (als)  ; 
f =cell2struct (answer (2 ) , ' d' , 1) ; 
a2s=f . d; 

a2=str2num (a2s)  ; 

OKf lag=TRUE; 
if  ( isempty (al ) ) , 

h=errordlg ( [als , '  is  not  a  number '],' Bum  Data  En¬ 
try'  )  ; 

OKflag=FALSE; 
waitf or (h) ; 
elseif  (isempty (a2 )} , 

h=errordlg ( [a2s , '  is  not  a  number '],' Bum  Data  En¬ 
try'  )  ; 

OKflag=FALSE; 
waitf or (h) ; 

end 

end 
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Appendix  A 


if  a2==0 / a2=l ; end; 
thru„attenuation=al  +  20  *loglO {a2 ) 


%  ask  user  name  of  file  for  direct: 
pname=f ilpre; 

[ f ildirect , pname]  =  uigetf ile ( [pname, asc Select  DIRECT 
file' ) 

if  ( ( f ilthru (1 ) ~=0)  &  ( fil thru ( length ( fil thru) )~=0) } , 
path_direct=pname ; 
else 

path_direct=f ilpre ; 

end 

%  enter  attenuation  for  DIRECT  file: 

OKflag^FALSE; 
while  OKflag==FALSE, 

prompt  =  {'Enter  dB  attenuation Enter  ratio  attenua¬ 
tion:  '  }  ; 

titlel  =  'Input  Attenuation  for  DIRECT  file'; 
lines =  1; 

def  =  { '20' , '!■ } ; 

answer  =  inputdlg (prompt , ti tlel , lines , def ) 

f  =-cell2struct  (answer  ( 1 )  ,  ' d '  ,  1 )  ; 
als=f . d; 

al  =  str2num (als )  ; 

f =cel 12s true t (answer (2 )  ,  ' d '  ,  1 )  ; 

a2s=f.d; 

a2=str2num (a2s )  ; 

OKflag=TRUE; 
i f  (isempty (al ) ) , 

h=errordlg ( [als , '  is  not  a  number '],' Bum  Data  En¬ 
try'  )  ; 

OKflag=FALSE; 
wait for (h) ; 
elseif  ( isempty (a2 )) , 

h=errordlg ( [a2s , '  is  not  a  number Bum  Data  En¬ 
try'  )  ; 

OKf lag=FALSE; 
waitf or (h) ; 

end 

end 

if  a2 = = 0 , a2 =1 ; end ; 

direct_.attenuation=al+20*logl0  (a2 ) 


%  Read  in  files: 

filname  =  streat ( f ilpre, deblank ( filthru) ) ; 
[x,y]  =  ascread (filname) ; 

scalefactor_thru=10''  (thru_attenuation/20)  ; 
y=y  *  (10^ ( thru_attenuation/20 ) ) ; 
filname  =  streat ( f ilpre, deblank ( fildirect) ) ; 
[xb,yb]  =  ascread (filname) ; 
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yb=yb  *  (10'"  (direct__attenuation/20  )  )  ; 

%  Check  if  sample  rates  are  the  same: 
if  abs { (x (2) -xb (2 ) ) /x (2 } )  >  9e-7, 

h=warndlg ([' Sampling  rates  are  different.  Press  any  key 
to  interpolate  time  domain  waveforms N / 'Warning' ) ; 
OKflag=FALSE; 
waitfor (h) ; 

%  interpolate  time  domain  waveforms: 

%  assume  thru  is  slower  rate  than  direct: 

%  pad  direct  (yb)  out  ot  end  of  thru  (y) .  Keep  sample 
rate : 

tstop=x (length (x) ) ; 
dt=xb(2 ) -xb (1); 
k=length (xb) ; 
while  xb (k) +dt<=tstop, 
k=k+l ; 

xb (k) -xb (k -1 ) +dt ; 
y b  ( k )  =  0  ; 

end 

%  interpolate  y  out  to  same  M  pts  of  yb  (direct) : 
method- ' spline' ; 
y=interpl (x,y,xb, method) ; 
x=xb; 


%  trunkate  to  2^14=16384: 
x=x(l:2'"14)  ; 
xb=xb(l:2"14) ; 
y=y(l:2-14) ; 
yb=yb(l:2'"14)  ; 


end 

disp ( ' Interpolation  done ' ) ; 


%  take  derivative  of  direct  pulse: 
dt=xb{2)-xb(l) ; 
for  i=2 ; length (yb) , 

dy (i) = (yb{ i) -yb (i-1 ) ) ;  %  dont  divide  by  /dt; 
end 

dy(l)=0; 

[x,y] =baseline2 (x,y, 10 , 138) ;  %  do  zero  adjust 
%  scale  direct  deriv: 

y_dir_deriv  =  dy  *  max(abs(yb))  /  max (abs (dy) ) ; 
f a=f f t (y ) ; 

fb=f ft (y_dir„deriv) ; 
f c=fa . /fb; 
yc=if f t ( fc ) ; 
npts2=arraylen ( fc) ; 
b=f ix (npts2/2 ) ; 
npts=b+l ; 
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%  do  a  smooth - 


f c_smooth-f c ;  %  er,  dont  do  a  smooth 
% [ fc_smooth] =smooth_complex ( fc, fix ( . 005*length { fc) } ) 

[I , J] =size (x) ; 
a=x(2) -x(l) ; 
fx=(0:npts-l) / ( {J*a) ) ; 
fx2= (0:npts2-l) / ( (J*a) ) ; 

fbD=fft(yb);  %  take  fft  of  direct 

%  plot  out  direct  &  fft,  direct-deriv  &  fft: 

Idir=l; 

displayf req=ll . 9e9 ; 

displayfreq=min (displayfreq, fx (length ( fx) ) ) ; 
while  fx (Idir) <displayfreq,  Idir=Idir+l;  end 
Idir=Idir-l ; 


[pk, fw, risl , xjunk,yjunk]  =  stats (x, abs (y) ) ; 

[pk, fw, ris2 , xjunk,yjunk]  =  stats (xb, abs (yb) ) ; 
figure, supti tie  (['Direct  f ile: ' , f ildirect] ) ; 
subplot (2,1,1) ,plot (xb,yb) ;grid  on; 
xlabe3('time  (S)');  ylabel ( 'Voltage  (V) ' ) ; 

%ti tie (' Response  of  PSPL  PulserM ; 

%legend ( ' thru ' , " direct ' ) ; 
v=axis ; 

a=v(l)  -I  (v(2)-v(l}  )*.55; 

bb  -  v(3)  +  ( (v(4)-v(3) )  *  0.95); 

text (a, bb, [ ' rise  time:  ' , num2str ( risl /le-12 ) , '  pS']) 
%a=v(l)  -t  (v(2)-v(l)  )  *.55; 

%b=v(3)  +  (v(4)-v(3) ) *.90; 

%  text (a, b,  [' fall  time(green):  ' , num2str (risb/le-12 ) 
pS '  ]  )  ; 

subplot (2,1,2), plot ( fx (1 : Idir) , 

20*logl0 (abs (fbD(l: Idir) ) ) ) ;  grid  on; 

xlabel (' frequency  (Hz)');  ylabel('(V/Hz)'); 

% [xf , yf ] =f f t_ps (x,y) ; 

% [xbf ,ybf ] =f f t_ps (xb,yb) ; 

%xcf =xf ; 

figure, suptitle  (['Direct-deriv  : ' , f ildirect ] ) ; 
subplot (2,1,1) ,plot (xb,y„dir_deriv) ;grid  on; 
xlabel ('time  (S) ' ) ;  ylabel ( 'normalized  (V/S) ' ) ; 
%title ( 'Response  of  PSPL  Pulser'); 

%legend ( ' thru ' , ' direct ' ) ; 
v=axis ; 

a=v(l)  +  (v(2)-v(l) )*.55; 
bb=v(3)  +  (v(4)-v(3) ) *.95; 

text (a, bb, [ ' rise  time:  ', num2str {ris2 /le-12 ), '  pS'3) 
%a-v(l)  +  (v(2) -v(l) ) *.55; 

%b=v(3)  +  (v(4) -v(3) ) *.90; 

%  text (a, b, [' fall  time{green):  ', num2str (risb/le-12 ) 
pS']  )  ; 
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subplot (2,1,2) , plot ( f X ( 1 : Idir ) , 

20*logl0 (abs ( fb (1 : Idir) ) ) ) ;  grid  on; 

xlabel  (' frequency  (Hz)');  ylabel(MV/Hz)'); 
%  [xf  ,yf  ]  =  f  f  t_ps  (x,y)  ; 

% [ xbf , ybf ] =  f  f  t_ps ( xb , yb ) ; 

%xcf =xf ; 


% 

%  let  user  do  filtering  of  the  transfer  function:- 


%[xl,yl]  = 

ginl ( fx, 20*logl0 (abs ( f c_smooth ( 1 :b+l ) ) ) , ' Hz ' , ' dB' ) 
%y2=10^ (yl/20) ; 

% [xt ,yt] =f loor_f irst (fx, abs ( f c ( 1 :b+l ) ) , xl ,y2 ) ; 

%%  take  the  max  of  y2  and  fc(i) : 

%for  i=l:k-l, 

%  if  y2>abs ( fc (i) ) 

%  fc„new(i)  =  (y2  ./  abs{fc(i)))  .*  fc(i); 

%  end 
%end 

%%  now  do  back  half: 

%for  i=m-(k“l)  :  m, 

%  if  y2>abs ( fc (i) ) 

%  fc_new(i)  =  (y2  . /  abs(fc(i)))  .*  fc(i); 

%  end 
%end 


[pk, fw, ris ,x,yjunk]  =  stats (x, abs (y) ) ; 

[pk, fw, risb, xb,yjunk]  =  stats(xb,abs(yb)); 
figure , suptitle  (filthru); 
subplot (2,2,1) ,plotyy (x,y,xb,yb) ;grid  on; 
xlabel ('time  (S)');  ylabel ( 'Voltage  (V)'); 
title ( 'Response  of  PSPL  Pulser'); 

%legend ( ' thru ' , ' direct ' ) ; 

%v=axis ; 

%a=v(l)  +  (v(2)-v(l) ) *.55; 

%b-v(3)  +  (v(4) -v(3) ) *.95; 

%  text (a, b, [' fall  time  (blue):  ' , num2str (ris/ le"12 ) , ' 
pS'] )  ; 

%a=v(l)  +  (v(2) -v(l) ) * . 55; 

%b=v(3)  +  (v(4)-v(3) ) *.90; 

%  text (a, b, [' fall  time (green):  ' , num2str (risb/le-12 ) , ' 
pS'  ]  )  ; 

%subplot (2 , 1 , 1) ,plot (x,y) ;  grid  on; 

%[xf ,yf]-fft_ps(x,y) ; 

% [xbf , ybf] =f f t_ps (xb,yb) ; 

%xcf  =:xf  ; 


f2=fa(l:b+l) ; 
f2 (l)=f2  (1) /2; 
f3=fb(l:b+l) ; 
f3 (l)=f3  (1) /2; 
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J=810;  %  J=200; 

sulqplDt(2,2,2)  ,Ed£t{fx(l:J)  ,20*IjGgl0{al:B(f2(l:J) ) )  ,fx(l:J)  ,20*lcgl0(ate(f3  {1:J) ) } )  ;grid 
on; 

xlabel  (' frequency  (Hz)M;  ylabel{'(dB)'); 
title ('FFTs  of  PSPL  data'); 
legend ( ' thru ' , ' direct '  ,  0 )  ; 

%subplot (2,2,3) ,plot(fx,20*logl0(abs(fc(l:b+l) ) ) ) ;grid  on; 
%xlabel (' frequency  (Hz)');  ylabel ( ' (dB) ' ) ; 

%title('FFT  of  PSPL  data'); 

%  plot  smoothed  curve: 

sukpaot(2,2,3)  ,plot(fx,20*logl0(abs(fc(l:bfl) ) )  ,fx,20*logl0(^(fc__srrnDfh{l:bfl) ) ) ) 

xlabel (' frequency  (Hz)');  ylabel ( ' (dB) ' ) ; grid  on; 

title  ( 'Transfer  funct'); 

legend ( ' original ' , ' smoothed ' , 0 ) ; 

%  plot  filtered  curve  but  only  up  to  10  GHz: 

1=1 ; 

while  fx(I)<10e9, 

1=1+1; 

end 

1=1-1; 

fc_new  =  fc_.smooth; 

SLiplot(2,2,4)  ,plot(f>c(l:I)  ,20*]c^0(^(fc_3TOOdi(l:I) ) )  ,fec(l:I)  ,20*logl0(dx(fcnew(l:I) ) ) ) 
xlabel (' frequency  (Hz)');  ylabel { ' (dB) ' ) ; grid  on; 
title ( 'Transfer  funct'); 
legend ( ' smoothed ' , ' truncated ',0); 

f 2 = f c_new ( 1 : b+ 1 ) ; 
f2  (l)=f2  (1) /2; 

f igure;plot ( fx2 , 20*logl0 (abs ( fc_new) ) ) ;grid  on; 

xlabel (' frequency  (Hz)');  ylabel {' (dB) ') ; 

title ( 'Transfer  Function  of  thru  using  PSPL  data'); 

yc_new= i f  f  t ( f  c_new ) ; 

figure; plot (x, real (yc_new) ) ; grid  on; 

xlabel ('time  (S)');  ylabel ( 'normalized  response  O'); 
title ( 'Impulse  Response  of  thru  using  PSPL  data'); 

%  Begin  the  plots  that  we  want: 

figure, %title  (filthru); 
plotyy (x, y , xb, y_dir_deriv) ;grid  on; 
xlabel ('time  (S) ' ) ;  ylabel ( 'Voltage  (V)  '  )  ; 
title(['PSPL  Pulser',',  ',filthru]); 

displayFreq=19 . 9e9 ;  c=l; 

while  (  (c<=npts)  &  (fx (c) <displayFreq)  ), 
c=c+l; 

end 
c=c-l ; 

f2=fa(l:b+l) ; 
f2(l)=f2(l)/2; 
f3=fb(l:b+l) ; 
f3 (l)=f3 (1) /2; 
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f 4=f c_smooth ( 1 :b+l ) ; 
f4 {l)=f4  (1) /2; 


J=810;  %  J=200; 
figure; 

plot(fx(l:c) .20*logl0(abs(f2 (l:c) ) ) ,fx(l:c) ,20*logl0 (afcs{f3 (1 :c) ) ) ) ;grid 
on; 

xlabel ( ' frequency  (Hz ) ' ) ;  ylabel ( ' (dB) ' ) ; 
title (['FFTs  of  PSPL  data  : ' , f ilthru] ) ; 
legend ( ' thru ' , ' direct ' , 0 ) ; 

displayFreq=ll . 9e9 ;  c=l; 

while  {  (c<=npts)  &  ( fx (c ) <displayFreq)  ), 
c=c+l; 

end 
c^c  1; 


figure; plot (fx (1 : c) , 20*logl0 (abs {f4 (1 ; c) ) ) } ;grid  on; 
xlabel (' frequency  (Hz)');  ylabel('(dB)'); 
title  ([ 'Transfer  function  of  , f ilthru] ) ; 


% 

% 

[newf ile, newpath]  = 

u  tf  ile  {  [path__thru ,  '  animini  t  ,  xf  r  '  ]  ,  '  Save  file  name' )  ; 

if  {  (newf  ile  ( 1 ) -=0 )  &  (newf  i-1  e  ( length  (newf  ile)  )  ~=0 ))  , 
ascwru.te„corrplex(fx2, fcjiew, filthru,  'legend' ,  'Hz' ,  'au' ,  [newpath, neArf:ilel )  ; 
end 


%  Compensateit3 . m 
%  4/99  N  Tesny 

%  corrects  multi-rate  pulser  data 

%  reads  in  already-computed  transfer  function  from  file 

%  Set  variables; 

FALSE=  0 ; TRUE=~FALSE ; 

%yes=true;no=f alse; 

filpre  =  'g:\uwb\other_subjects\lpDecon\'; 

%filpre  =  'c:/dascw/' 

%ItemNumber  =  6; 

%  files  to  calculate  transfer  function: 
%fildirect='drc2pspl.asc' ; 

%f ilthru= 'pspllpLogP. asc ' ;  %  Ip 
%f ildirect= ' drc2pspl . asc ' ; 

%f ilthru= 'psplT220SS . asc ' ;  %  t2 
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%  ask  user  name  of  file  for  thru: 
pname-filpre; 

[  f  ilthru ,  pname]  uigetf  ile  (  [pname asc  ^  ^  Select  THRU 
file" ) 

if  ( {filthru(l)~=0)  &  { fil thru ( length ( fil thru) )-=0) ) , 
path__thru=pname  ; 
else 

path_thru=f ilpre; 

end 

%  enter  attenuation  for  THRU  file: 

OKflag=FALSE; 
while  OKf lag==FALSE, 

prompt  =  {'Enter  dB  attenuation: Enter  ratio  attenua¬ 
tion:  ' }  ; 

titlel  =  'Input  Attenuation  for  THRU  file'; 
lines-  1; 

def  =  {'O', '2'}; 

answer  =  inputdlg (prompt , titlel , lines , def ) 

f~ cel 12 struct (answer ( 1 } , ' d' , 1 ) ; 
als=f .d; 

al=str2num (als )  ; 
f =cell2struct (answer (2 ) , 'd' , 1) ; 
a2s=f . d; 

a2=str2num (a2s) ; 

OKflag=TRUE; 
if  (isempty (al ) )  , 

h=errordlg ( [als , '  is  not  a  number '],' Bum  Data  Entry'); 
OKflag=FALSE; 
wai tfor (h) ; 
elseif  ( isempty (a2 )) , 

h=errordlg( [a2s, '  is  not  a  number' ] , 'Bum  Data  Entry' ) ; 

OKflag=FALSE; 

wait for (h) ; 

end 

end 

if  a2==0,a2=l;end; 
thru„attenuation=al+20  *loglO (a2 ) 


%  ask  user  name  of  file  for  direct: 
pname=f ilpre; 

[ f ildirect , pname]  =  uigetf ile ( [pname asc '],' Select  DIRECT 
file' ) 

if  ( ( f ildirect (1)“=0)  &  ( f ildirect (length ( f ildirect ) )~=0) ) , 
path_direct=pname ; 
else 

path„direct=f ilpre; 

end 

%  enter  attenuation  for  DIRECT  file: 

OKflag=FALSE; 
while  OKflag==FALSE, 

prompt  =  {'Enter  dB  attenuation: Enter  ratio  attenua¬ 
tion:  '  }  ; 

titlel  =  'Input  Attenuation  for  DIRECT  file'; 
lines=  1; 

def  =  {  '20'  ,  '2  •  }  ; 
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answer  =  inputdlg {prompt, titlel, lines, def) 

f=cel 12s tract (answer (1) , "d^ , 1 ) ; 
als=f .d; 

al=str2num(als) ; 

f=cel 12s tract (answer (2) , 'd' , 1 ) ; 

a2s=f . d; 

a2  =  str2nam (a2s)  ; 

OKflag=TRUE; 
if  (isempty (al) ) , 

h=errordlg ( [als , '  is  not  a  namber ' ] , ' Bam  Data  EntryM 
OKflag=FALSE; 
waitfor (h) ; 
elseif  ( isempty (a2 )) , 

h=errordlg ( [a2s , '  is  not  a  namber ' ] / ' Bam  Data  EntryM 

OKflag=FALSE; 

waitfor (h) ; 

end 

end 

if  a2-=0,a2=l;end; 
direct_attenaation=al+20  *loglO (a2 ) 

%  ask  aser  name  of  file  for  TRANSFER  FUNCTION: 
pname=^filpre; 

[ f i ] transfanc , pname]  =  aigetf ile ( [pname, .xfr Select 
file  to  ase  for  TRANSFER  FUNCTION') 
if  (( fil transfanc (1) ~=0 )  & 

( fil transfanc (length ( fil transfanc) ) ) ) 
path_transfanc=pname; 
else 

path_.  transfunc  =  f  ilpre; 

end 


%  Read  in  files: 

fjJname  =  strcat (path_thra, f ilthra) ; 

[x,y]  =  ascread { f ilname) ; 

y  y  *  10^  { thru_.attenaation/20 )  ;  %  scale  for  attenuation 
%  plot  thru: 

[pk, fw,  ris, x,yjank]  =  stats(x,abs(y)); 
figure; 

plot (x,y) ; grid  on; 
t i t] e  ( 'Captured  Signal'); 
v-axis ; 

a-v(l)  +  (v(2)-v(l) ) *.55; 
b-v(3)  +  (v(4)-v{3))*.95; 

text (a, b, [' rise  time  (blue):  ' , num2str (ris/le- 12 ) , '  pS']); 
%  read  in  transfer  function: 

filnamt  =  strcat (path_transfunc , filtransfunc) ; 

[x„tf ,y„real ,y_imag]  =  ascread_complex ( filnamt ) ; 
y„tf =y„real+ j  *y_imag; 

%  plot  tf: 

figure; plot (x„tf , 20*logl0 (abs (y_tf ) ) ) ; grid  on; 
title ( 'Transfer  Function'); 


%  Read  in  direct  file: 
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f  ilnameDirect  =  strcat  (path..direct ,  f  ildirect )  ; 

[x_dir ,y_dir]  =  ascread ( f ilnameDirect ) ; 

y_dir  =  y_dir  *  lO'"  (direct„attenuation/20)  ;  %  scale  for 
attenuation 

y_dir„f f t=f f t (y_dir ) ; 


npts2__d=arraylen  (y_dir_f  f  t )  ; 
b=f ix {npts2_d/2 ) ; 
npts_d=b+l; 

[I, J] =size (x_dir ) ; 
a=x_dir ( 2 ) -x_dir (1); 

fx_d= (0 :npts_d-l) / ( (J*a) } ;  %  half  freq  data 
%  fx2_d=(0:npts2„d-l)/ ( (J*a) ) ;  %  full  freq  data 


%  plot  direct  and  its  fft: 
figure, %suptitle ( filthru) ; 
subplot (2,1,1) , 
plot (x_dir ,y_dir} ; 

grid  on;  xlabel('time  (S) ' ) ;  ylabel { ' Vol tage  (V) ' ) ; 
title (' Response  of  Direct'); 


Idir=l; 

displayf req=10e9 ; 

displayfreq=min (displayfreq, fx„d ( length ( fx_d) ) ) ; 
while  fx_d (Idir ) <displayfreq,  Idir=Idir+l;  end 
Idir=Idir-l ; 

subplot (2,1,2)  , 

plot ( fx_d (1 : Idir) , 20*logl0 (abs (y_dir_f f t (1 : Idir) ) ) ) 
grid  on;  xlabel (' frequency  (Hz)');  ylabel (' (dB) ') ; 
title  ( 'FFTs' ) ; 


%  Check  to  see  if  we  need  to  interpolate: 

%  ie,  if  sampling  rates  are  different: 

%  take  fft: 
fa=f f t (y) ; 
npts2=arraylen ( fa) ; 
b=f ix (npts2/2 ) ; 
npts=b+l ; 

[I, J] =size (x) ; 
a:=x(2)  -x(l)  ; 

fx= (0 :npts-l) / ( (J*a) ) ;  %  half  freq  data 
fx2= (0:npts2-l) / ( (J*a) ) ;  %  full  freq  data 

if  abs( (x_tf (2)-fx2(2) ) )/fx2(2)  >  9e-7, 

h=warndlg ([' Sampling  rates  are  different.  Press  OK  to 
interpolate  the  transfer  function' ], 'Note' ) ; 
waitfor (h) ; 
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%  break  the  fft  and  tf  in  half: 
npts_full  =  length (fa); 
npts_half  =  ceil (npts_full/2+ . 5 ) ; 
fa_half=fa (1 :npts_half ) ; 

fa_half (1 ) =fa_half { 1 ) /2 ;  %  divide  dc  value  in  half 

npts_full„tf  =  length (y„tf ) ; 
npts_half__tf  =  ceil  (npts_full_tf /2+ .  5 )  ; 
tf„half  =  y_tf (1 :npts_half_tf ) ; 

tf_half(l)  =  tf_half (1) /2;  %  divide  dc  value  in  half 
x„tf_half  -  x_.tf  ( 1 : npts_half_tf )  ; 

%  do  interpolation  of  tf: 

%method= ' linear' ; 
method= ' nearest ' ; 

y..  tf_interp  =  interpHx_tf__half ,  tf_half;  fx,  method); 

%  now  reconstruct  inturpolated  wf  to  full  length: 

%  remember  that  mirrored  side  is  complex  conjugate! 

yTf  InterpFull=ones  (1  ,  npts_.full )  ; 

yTf  InterpFull  (1 : npts._hal f )  =y_tf„interp; 

for  i=2 : npts„half -1 , 

yTf  InterpFuI  1  (npts...  full +2-i )  •  conj  (yTf  InterpFull  (i  )  ) 

end 

yTf InterpFuI 1 ( 1 )  =  yTf InterpFull ( 1 )  *  2;  %  mult  dc  value 
by  2 
el  se 

npts._.full  -•  length  (fa); 
npts„half  =  ceil(npts._full/2+.5); 
yTf InterpFull  =  y_tf; 

end 

%  Now  apply  transfer  function  to  fft: 
fb=  fa . /yTf InterpFull ; 

%  Filter  it  out: 
fi 1 terflag-FALSE; 
fpHi=3e9 ; 
f sHi =6e9 ; 
fpT.o--200e6 ; 
f sLo=50e6 ; 

delHi=-80 ; delLo=- 20 ; 
fpDi^l0e9; 

fsDj=20e9;  delDi--40; 

%  begin  filtering 

loop _ 

while  f ilterf lag==FALSE, 

OKflag=FALSE; 
while  OKf lag— FALSE, 

prompt  =  {'Enter  F-passband  in  GHz:', 'Enter  F- 
stopband  in  GHz: ' , 'Enter  stopband  loss  in  dB  (delta) ; ' } ; 
titlel  =  'HIGH  Frequency  Filtering'; 
lines=  1; 

def  =  {num2str ( fpHi/le9 ) , num2str ( f sHi/ 

le9 ) , num2str (delHi) } ; 

answer  =  inputdlg (prompt , titlel , lines , def ) 

f= cel 12 struct (answer ( 1 ) , ' d' , 1 ) ; 
als=f .d; 
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fpHi=str2num (als) ; 
f =cell2struct (answer (2 ) , 'd' , 1) ; 
a2s=f . d; 

f sHi=str2num (a2s ) ; 

f=cell 2 struct (answer (3 ) , ^  d ' . 1) ; 

a3s=f . d; 

delHi=str2num (a3s ) ; 

OKflag=TRUE; 
if  (isempty ( fpHi) ) , 

h=errordlg ( [als, '  is  not  a  number ^ Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  ( isempty ( f sHi ) ) , 

h=errordlg ( [a2s, '  is  not  a  number Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  (isempty (delHi) ) , 

h=errordlg ( [a3s , '  is  not  a  number Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

end 

end 

fpHi=fpHi*le9; 
f sHi=f sHi *le9 ; 

%  Let  user  enter  low  freq  numbers: 

OKflag=FALSE; 
while  OKf lag==FALSE, 

prompt  =  {'Enter  F-stopband  in  MHz:', 'Enter  F-pass- 
band  in  MHz:', 'Enter  stopband  loss  in  dB  (delta);'}; 
titlel  =  'LOW  Frequency  Filtering'; 
lines=  1; 

def  =  {num2str ( f sLo/le6 ) , num2str { fpLo/ 

le6 ) , num2str (delLo ) } ; 

answer  =  inputdlg (prompt , titlel , lines, def ) 

f =cel 12s true t (answer ( 1 ) , ' d ' , 1 )  ; 
als=f . d; 

f sLo=str2num (als) ; 

f=cel 12 struct (answer (2 ) , ' d' , 1)  ; 

a2s=f . d; 

fpLo=s tr2num (a2s ) ; 
f =cell2struct (answer (3 ) , 'd' , 1)  ; 
a3s=f . d; 

delLo=str2num (a3s) ; 


OKflag=TRUE; 
if  (isempty (fpHi) ) , 

h=errordlg ( [als, '  is  not  a  number '],' Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  ( isempty ( fsHi )) , 
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h=errordlg { [a2s , '  is  not  a  number M Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  (isempty (delHi) ) , 

h=errordlg ( f a3s , '  is  not  a  number Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

end 

end 

fpLo=fpLo*le6 ; 
f sLo=f sLo*le6 ; 

%  let  user  enter  limits  for  direct,  waveform  filter: 

OKflag=FALSE; 

while  OKf lag==FALSE, 

prompt  =  {'Enter  F-passband  in  GHz:', 'Enter  F-* 
stopband  in  GHz:', 'Enter  stopband  loss  in  dB  (delta);'}; 

titlel  =  'HIGH  Frequency  Filtering  For  DIRECT  Wave¬ 
form'  ; 

lines=  1; 

def  =  {num2str ( fpDi/le9 ) , num2str ( f sDi/ 

le9 ) , num2str (delDi) } ; 

answer  =  inputdlg (prompt , ti t lei , 1 i nes , def ) 

f=cel 12s true t (answer ( 1 ) , ' d ' , 1 ) ; 
als=f . d; 

fpDi-str2num (als ) ; 
f =cell2struct (answer (2 ) , ' d ' , 1 ) ; 
a2s=f  .  d; 

fsDi=str2num (a2s ) ; 
f-cell2struct (answer (3 ) , ' d' , 1 ) ; 
a3s=f . d; 

delDi=str2num (a3s ) ; 

OKflag=TRUE; 

if  ( isempty ( fpHi ) ) , 

h=errordlg ( [als, '  is  not  a  number '],' Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  ( isempty ( fsHi )) , 

h=errordlg ( [a2s , '  is  not  a  number '],' Bum  Data 

Entry' ) ; 

OKflag=FALSE; 
waitfor (h) ; 

elseif  (isempty (delHi )) , 

h=errordlg ( [a3s , '  is  not  a  number '],' Bum  Data 

Entry' ) ; 

OKflag^FALSE; 
waitf or (h) ; 

end 

end 

fpDi=fpDi *le9 ; 
f sDi=f sDi *le9 ; 


%  Generate  filter  from  user  specs: 
for  i=l : length ( fx) , 
if  fx(i)<=fsLo, 
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filtl (i)=delLo; 

elseif  fsLo<fx(i)  &  fx(i)<fpLo, 

filtl  (i)  =  delLo  *  (fpLo-fx(i))  /  ( fpLo- f sLo) ; 

elseif  fpLo<=fx(i)  &  fx(i)<=fpHi, 
filtl(i)=0; 

elseif  fpHi<fx{i)  &  fx(i)<fsHi, 

filtl  (i)  =  delHi  *  (fpHi-fx(i))  /  (fpHi-fsHi)  ; 

elseif  fx(i)>=fsHi, 
filtl(i)  =  delHi; 

end 

filtl  (i)  =  10  ^  (f iltl (i) /20) ;  %  convert  from  dB  to 

number 

end 


%  now  reconstruct  filter  to  full  length: 
f ilt_full=zeros (l,npts_full) ; 
filt_full  (1  :npts_half )  ==filtl; 
for  i=2 : npts_half -1 , 

f ilt„full (npts_full+2“i)  =  f ilt_full (i ) ; 

end 

figure;plot (20*logl0 {filt_full) ) ; 

%  Scale  with  filter: 
fbFilt  =  f b  . *  filt_full; 

%  Do  if ft: 

y_new  =  if ft ( fbFilt ) ; 

%  check  waveform  for  complex  data: 
n=min  ( 40 ,  length  (y__new)  )  ; 
warningFlag^FALSE ; 
for  i=l:n, 

if  abs (imag (y_new (i) ) )  >  9e-7, 
warningFlag=TRUE ; 

end 

end 

if  (warningFlag==TRUE) , 

h=warndlg ([ 'Complex  data  is  present  in  the  time  domain 
waveform! ' ] , 'Warning' ) ; 
wait for (h) ; 

end 

figure ; plot (x, real {y_new) ) ; 

%  check  for  sample  rates  being  different: 
if  abs { (x_dir (2 ) ”X (2 ) ) ) /x (2 )  >  9e-7 , 

h=warndlg ([ 'Sample  rate  of  Direct  is  different  from 
Thru.  Press  OK  to  continue '],' Note ') ; 
waitf or (h) ; 

end 

%  Take  derivative  of  direct: 
dt=x_dir (2 ) -x_dir ( 1) ;  dy{l)=0; 
for  i=2 : length {y_dir) , 

dy(i)= (y_dir(i)-y_dir (i"l) ) ; 

end 

%  scale  direct  deriv: 

y_dir_deriv  =  dy  *  max (abs (y_dir ) )  /  max (abs (dy)); 
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dirDerivFFT  =  fft(y_dir_deriv); 

npts2_d=arraylen (dirDerivFFT) ; 
b=f ix (npts2_d/2 ) ; 
npts_d-b+l; 

[I, J] =size (x_dir)  ; 
a=x_dir ( 2 ) -x_dir (1); 

fx_d= (0:npts„d-l} / ( (J*a) ) ;  %  half  freq  data 
fx2_d= (0 : npts2_d-l ) / ( ( J*a} ) ;  %  full  freq  data 


%  Generate  filter  for  direct  from  user  specs: 
for  i=l : length ( fx_d) , 
if  fx_d (i } <=fpDi , 
filtD(i)  =  0; 

elseif  fpDi<fx„d{i)  &  fx_d (i) <fsDi , 

filtD(i)  =  delDi  *  ( fpDi-fx_d ( i ) )  /  (fpDi-fsDi); 
elseif  fx_d ( i ) >=f sDi , 
filtD(i)  =  delDi; 

end 

filtD{i)  =  10  ^  { f il tD (i ) /20 ) ;  %  convert  from  dB  to 

number 

end 


%  now  reconstruct  filter  to  full  length: 
f il t_fullD= zeros ( 1 , npts2_d) ; 
f ilt_fullD {1 :npts_d) =f iltD; 
for  i=2 : npts_d-l , 

f  ilt...fullD(npts2  .d4  2-i  )  ^  f  i  lt_f  ullD  (  i  )  ; 

end 


%  Sca]e  with  filter: 

di  rDerivFFT„.f  ilt  -  dirDerivFB'T  .*  filt„fullD; 
%  Do  if ft: 

y„dir_deriv„new  =  if f t (dirDerivFFT_f ilt ) ; 

%f igure;plot (x, real (y_dir_new) ) ; 

1  igure;plot  (x_dir  ,y__dir , 
x.  dir , real (y_dir_deriv_new) ) ;grid  on; 
title ( 'Filtered  direct  waveform'}; 
xlabel ( ' time  (S) ' ) ;ylabel ( 'V' ) ; 
legend ( ' direct ' , ' direct-deriv' ) ; 


%  Plot  corrected  and  direct: 

figure; plot (x, y_new, x_dir , y„dir_deriv_new) ; 

grid  on; xlabel (' time  (S) ' } ;ylabel ( 'V' ) ; 

title ( 'Results ' ) ; legend ( 'corrected' , ' direct  deriv' ) ; 

figure; plot ( fx,  20*logl0 (abs { fbFilt ( 1 : npts_half ))),... 

fx„d,  20*logl0 (abs {dirDerivFFT_f ilt (1 :npts_d) ) ) ) ; 
grid  on; xlabel (' freq  (Hz ) ' ) ; ylabel ( ' V/Hz ' ) ; 
title ( 'Results ' ) ; legend ( 'corrected' , ' direct  deriv' ) ; 


displayFreq=14 . 9e9 ;  b=l; 
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while  (  (b<=npts_half )  &  ( fx (b) <displ ayFreq)  ), 
b=b+l; 

end 

b=b-l; 

displayFreq=14 . 9e9 ;  bd=l; 

while  (  (bd<=npts_d)  &  ( f x_d (bd) <displayFreq)  ). 
bd=bd+l; 

end 

bd=bd-l; 

figure; plot ( fx (1 :b) , 20*logl0 (abs ( f iltl ( 1 : b) ) ) , ' r-' , . . . 
fx(l:b) ,20*logl0 (abs{fb(l:b) ) ) , 'k-' , . . . 
fx(l:b) , 20*logl0 (abs (fbFilt (l:b) ) ) , ^b-^ . . . 
fx_d(l:bd) , 20*logl0 (abs (dirDerivFFT_f ilt (l:bd) ) ) / 'g- 

^ )  ; 

grid  on;xlabel ( ' f req  (Hz) ' ) ;ylabel ( 'V/Hz ' ) ; 
title( 'Results' ) /legend ( 'filter' , 'unfiltered' , 'filtered' , 'direct 
deriv' ) ; 

filterflag=TRUE; 

button  =  questdlg ( ' Are  filter  limits  OK?',.,. 

'Continue  Operation' , 'Yes ' , 'No' , 'Help' , 'No' ) 
if  strcmp (button, ' Yes ' ) 
disp ( ' OK' ) 

elseif  strcmp (button, 'No' ) 

disp ('Redo  filter  operation') 
f ilterf lag=FALSE; 

close  4 /close  5 /close  6 /close  7 /close  8 /close  9; 
elseif  strcmp (button, ' Help' ) 

disp ('Sorry,  no  help  available') 

end 


end;  %  end  of  loop  for  filterflag 


[pkl , fwl , risl , xjunk,yjunk]  =  stats (x, real (y„new) ) ; 
rtl=num2str (round (risl/le-12 ) ) ; 

[pk2 , fw2 , ris2 , xjunk,yjunk]  =  stats (x,y); 
rt2=num2str (round (ris2/le-12 )  )  ; 

[pk3 , fw3 , ris3 , xjunk,yjunk]  = 

stats (x_dir, real (y_dir_deriv_new) ) ;  rt3=num2str (round (ris3 / 
le-12) ) ; 


%  normalize  and  plot  on  same  graph: 
figure, suptitle ( filthru) ; 
subplot (2,1,1) , 

[xAlignl,y_newAlign]  =align__peak {x,y_new)  ; 

[xAlign2  ,yAlign]  =align__peak  (x,y)  ; 

[xAlign3,dyAlign]  =align_peak  (x_dir  ,y_dir_deriv_new)  ; 


timbeg=-2e-9 ; timend=2e-9 ; 
if  strncmpi (filthru, 'av' , 2 ) , 

timbeg=-9.9e-9/timend=20e-9;  %if  avtek,  display  longer 
window 
end 

if  strncmpi ( filthru, 'hhlp' , 4 ) , 
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timbeg=-9.9e~9;  timend=20e--9;  %if  Ip,  display  longer 
window 
end 

timbeg=-20e-9; timend=20e-9 ; 

[xAlignl ,y_newAlign] =prepst (xAlignl ,y_newAlign, timbeg, timend) ; 
[xAlign2  ,yAlign]  =prepst  (xAlign2  ,  yAlign,  timbeg,  timend)  ; 
[xAlign3  ,  dyAlign]  =prepst  (xAlign3  ,  dyAlign,  timbeg,  timend)  ; 

if  abs {max (yAlign) ) <abs (min (yAlign) ) , 
yAlign=-yAlign; 

end 

plot (xAlignl ,y_newAlign/max (y„newAlign) , ^g' , xAlign2 , yAlign/ 
max (yAlign) , . ♦ . 

'b' ,xAlign3 , dyAlign/max (dyAlign) , 'r ' ) ; 
grid  on;  xlabel('time  (S)');  ylabel ( 'Vol tage  (V)  M ; 

%title  ( 'Response  of  PulserM; 

legend  ([ 'comp,  rt=',rtl,'  ps '  ]  ,  T 'meas ,  rt=',rt2,'  psM/f'd/ 
dt ,  rt= ' , rt3 , '  ps ' ] ) ; 

1=1; 

displayf reg=10e9 ; 

displayfreq=min (displayfreq, fx ( length ( fx) ) ) ; 
while  fx (I) <displayfreq,  1=1+1;  end 
I=I~1; 


Idir=l; 

displayf req=10e9 ; 

displayf  req=min  (displayf  req,  fx_d  ( length  { fx_,d)  )  )  ; 
while  fx„d (Idir ) <displayfreq,  Idir=Idir+l;  end 
Idir=Idir-l ; 


subplot (2,1,2) , 

plot(fx(l:I) ,20*logl0(abs(fa(l:I) )),... 
fx{l:I) ,20*logl0{abs(fbFilt(l:I) )),.., 

fx„d(l:Idir)  ,  20*logl0  (abs  (di  rDer.ivFFT_fil  t  (1 :  Idir )  )  )  )  ; 

grid  on;  xlabel (' frequency  (Hz) M ;  ylabel (' (dB) ') ; 
title  ( ' FFTs ' ) ; 

legend { 'meas ' , " comp' , ' d/dt ' , 0 ) ; 


Appendix  B.  Examples  of  Compensation 


This  appendix  shows  examples  of  the  original  data  and  the  intermediate 
results  as  the  data  are  processed  according  to  the  flow  charts  discussed  in 
figures  6  and  7.  Figures  B-1  through  B-4  show  the  direct  4015  signal  re¬ 
corded  on  the  scope  as  well  as  that  recorded  through  the  T3  antenna  pair. 
Figures  B-5  through  B-8  show  the  fast  Fourier  transforms  of  these  recorded 
signals.  Figures  B-9  through  B-12  show  the  calculated  transfer  functions  of 
the  4015  through  the  four  antermas,  and  figures  B-13  through  B-21  show  the 
reconstructed  signals.  The  reconstructions  are  plotted  over  the  original  sig¬ 
nal  for  comparison,  and  these  compare  well.  The  complete  matrix  of  sources 
used  during  the  experiments  is  shown  in  table  B-1. 


Table  B-1.  Antennas  and  sources  studied. 


Source 

Antennas 

T3 

T2 

LP 

EMCO 

PSPL4015C 

X 

X 

X 

X 

HH2 

X 

X 

X 

X 

Avtech 

X 

X 

X 

X 
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Figure  B-1 .  Signals 
from  PSPL  source  and 
T2  antenna  pair. 
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Figure  B-2.  Signals 
from  PSPL  source  and 
EMCO  antenna  pair. 
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Figure  B-3.  Signals 
from  PSPL  source  and 
log  periodic  antenna 
pair. 
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Figure  B-4.  Signals 
from  PSPL  source  and 
T3  antenna  pair. 
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Figure  B-5.  FFTs  of 
signals  from  PSPL 
source  and  T2  antenna 
pair. 
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Figure  B-6.  FFTs  of 
signals  from  PSPL 
source  and  EMCO 
antenna  pair. 
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Figure  B-7.  FFTs  of 
signals  from  PSPL 
source  and  log 
periodic  antenna  pair. 


Figure  B-8.  FFTs  of 
signals  from  PSPL 
source  and  T3  antenna 
pair. 
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Figure  B-9.  Transfer 
function  of  T2  antenna 
pair. 


Figure  B-10.  Transfer 
function  of  EMCO 
antenna  pair. 
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Figure  B-ll.  Transfer 
function  of  log 
periodic  antenna  pair. 


Figure  B-12.  Transfer 
function  of  T3  antenna 
pair. 
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Figure  B-13.  1 
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Figure  B-14. 
Reconstruction  of 
signal  from  Avtech 
source  through  LP 
antenna  pair. 
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Figure  B-17.  Jitter  in 
Avtech  source. 


Figure  B-18.  Multi-rate 
reconstruction  of 
signal  from  HH2 
source  through  T2 
antenna  pair. 
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Figure  B-19.  Multi-rate 
reconstruction  of 
signal  from  HH2 
source  through  T3 
antenna  pair  measured 
with  an  SCDIOOO 
oscilloscope. 
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Figure  B-20.  Multi¬ 
rate  reconstruction  of 
signal  from  HH2 
source  through 
EMCO  antenna  pair 
measured  with  an 
SCDIOOO 
oscilloscope. 
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